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Attention is drawn to the role played by the size of the system in the thermodynamic analysis 
of particle yields in relativistic heavy ion collisions at SIS energies. This manifests itself in the 
non-linear dependence of K + and K~ yields in AA collisions at 1 - 2 ^4-GeV on the number of 
participants. It is shown that this dependence can be quantitatively well described in terms of a 
thermal model with a canonical strangeness conservation. The measured particle multiplicity ratios 
(tt + /p, ir~/iv + , d/p, K + /tv + and K + /K~ but not 77/71"°) in central Au-Au and Ni-Ni collisions at 
O |i 0.8 - 2.0 A-GeV are also explained in the context of a thermal model with a common freeze-out 

temperature and chemical potential. Including the concept of collective flow a consistent picture 
of particle energy distributions is derived with the flow velocity being strongly impact-parameter 
dependent. 

I. INTRODUCTION 

> 

\ It was pointed out by Hagedorn |IJ some thirty years ago that thermal models overestimate the production of 
anti-He 3 in proton-proton collisions by seven orders of magnitude when the grand canonical ensemble is used in its 
0^ ' standard form |2|. The reason for this is that when the number of particles as well as the interaction volume are 
small one has to take into account the fact that the production of anti-He 3 must be accompanied by the production 
of another three nucleons with energy En in order to conserve the baryon number. Thus, the abundance will not be 
proportional to the standard Boltzmann factor given by 



but to 
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since three additional nucleons must be produced in order to conserve the baryon number. This suppresses the rate 
and introduces a cubic volume dependence. The original presentation of Hagedorn has been considerably developed 
and expanded in refs. [p|-|lT|. 

Recently it has become clear that a similar treatment should be followed for strangeness production in the GSI/SIS 
energy range fl2]| . This is not only due to the fact that the size of the system is small but mainly because the 
temperature is very low and particle numbers are small. The abundance of _ftT + -mesons is then given by 
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since the strangeness must be balanced either by an anti-kaon or by a hyperon. gi are the degeneracy factors and 
Ei the particle energies. This leads to a linear dependence of the K + density on the size of the system. Such a 
dependence has indeed been observed by the KaoS collaboration for if + -mesons [ p^JT^ I . 

In this paper we would like to explore this idea in detail. This volume dependence can now be tested for the first 
time by considering the data on impact parameter dependence which are now becoming available. 

In section II we review the thermal model with special emphasis on the canonical corrections due to the exact 
conservation of quantum numbers. In section III wc present a comparison with the experimental data from SIS. One 
part is devoted to a systematic study of central collisions of various systems, another part to a detailed investigation 
of the impact-parameter dependence. The last section is devoted to a discussion of our results. 
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II. CONCEPTS AND PREDICTIONS 



The exact treatment of quantum numbers in statistical mechanics has been well established for some time now [p| . 
It is in general obtained by projecting the partition function onto the desired values of the conserved charges by using 
the group theoretical methods (for a review see e.g. 0.) For our purpose we shall only consider the conservation 
laws related to the abelian U(l) symmetry group. In this case the formalism is particularly simple and leads to the 
following form of the canonical partition function for a fixed value of the conserved charge Q: 

where Z is obtained from the grand canonical (GC) partition function replacing the fugacity parameter Aq by the 
factor e lc ^, 

Z{T,V^) = Z GC (T,V,\ Q ^e^) (2.2) 

The particular form of the generating function Z in the above equation is model dependent. Having in mind the 
applications of the statistical description to particle production in heavy ion collisions we calculate Z in the ideal 
gas approximation, however, including all particles and resonances listed in [fl5|| . This is not an essential restriction, 
because, describing the freeze-out conditions, we are dealing with a dilute system where the interactions should not 
influence particle production anymore. We neglect any medium effects on particle properties. In general, however, 
already in the low-density limit, the modifications of resonance width or particle dispersion relation, in this particular 
for A and tt |l(|[l7|, cannot be excluded. For the sake of simplicity, we use classical statistics, i.e. we assume 
temperature and density regime so that all particles can be treated using Boltzmann statistics. 

In nucleus- nucleus collisions the absolute values of baryon number, electric charge and strangeness are fixed by 
the initial conditions. Modeling particle production in statistical thermodynamics would, in general, require the 
canonical formulation of all these quantum numbers. ^From the previous analysis p2|, however, it is clear that only 
strangeness should be treated exactly, whereas the conservation of baryonic and electric charges can be described by 
the appropriate chemical potentials in the grand canonical ensemble. 

Within the approximations des cribed above and neglecting the contributions from multi-strange baryons, the gen- 
erating function in equation (|2.l[) has the following form for a gas having zero total strangeness, S = 0: 

Zs(T,V,w,ii B ,<l>)=exp(N a= 4> + N a= ie i + + N a= -ie- i +) (2.3) 
where N s=Q ±i is defined as the sum over all particles and resonances having strangeness 0, ±1, 

N s =^±i = Y. Z l ( 2 - 4 ) 



and Z\ is the one-particle partition function defined as 

Z k = TTT m k T K 2{mk/T) cxp(b k fi B +qkVQ) (2.5) 

ZTT 

with the mass rrik, spin-isospin degeneracy factor the particle baryon number b k and electric charge q k - The 
volume of the system is V and the chemical potentials related with the charge and the baryon number are determined 
by fiQ and [Ib respectively. 

With the particular form of the generating function equations (2.3,2.4 2~5| ) the ^-integration in equation (2.1) can 



be done analytically giving the canonical partition function for a gas with total strangeness S in the following compact 
form 0: 

Zs(T, V, Mo) = z o(T, V, ^ b ,^q)I s {x) (2.6) 

where Zq — exp (Ns=o) is the partition function of all particles having zero strangeness and the argument of the 
Bessel function 



x = 2y/S 1 S- 1 . (2.7) 

with S±i = N s= ±i. 



2 



The calculation of the particle density nk in the canonical formulation is straightforward. It amounts to the 
replacement 

Zl i ► A fc Zl (2.8) 



of the corresponding one-pa rticl e partition function in equation (2.4) and taking the derivative of the canonical (C) 



partition function equation (2.1) with respect to the particle fugacity Afc 

(2.9) 



Q 

n k = A fc TTT-lnZQ(A fe ) 

OA k 



A fe = l 



As an example, we quote the result for the density of thermal kaons in the canonical formulation assuming the total 
strangeness of the system S = 0, 



Z X K Si h(x) 



(2.10) 



Comparing the above formula with the result for thermal kaons density in the grand canonical ensemble, n^P = 
(Z^ /V) exp (/ig /T), one can see that the canonical result can be obtained from the grand canonical one replacing the 
strangeness fugacity As = exp (fis/T) in the following way: 

nC K = nG K c ( x s^^^rn)- ( 2 - u ) 

\ y/SiS-! Iq{x) J 

In the thermodynamic limit both the canonical and the grand canonical formulation are equivalent. For a small 
system, however, the differences are large. This can be seen in the most transparent way when comparing two 
limiting situations: the large and small volume limit of equation (2.10) In the thermodynamic limit V — > oo the 
argument of the Bessel function x — > oo, thus 

km £M -> 1 (2.12) 

x^oa 1q[x) 

and the kaon density is independent of the volume of the system as expected in the grand canonical ensemble. On 
the other hand in the limit of a small volume we have 

L (x) x 

lim -^-4 -> - (2.13) 

and the particle density is linearly dependent on the volume. It is thus clear that the major difference between the 
canonical and the grand canonical treatment of the conservation laws appears through different volume dependence of 
strange particle densities. The relevant parameter, F$, which measures the deviations of particle multiplicities from 
their grand canonical result is determined by the ratio of the Bessel functions 

Fs = ^ (2.14) 
Io{x) 



with the argument x defined in equation (2.7). In Fig. 1 we show the canonical suppression factor Fs as a function of the 
argument x. To relate the initial volume of the system to the number of participants we use the approximate relation 
V ~ l.9irA part . The corresponding values of x at SIS, AGS and SPS energies are calculated with the baryochemical 
potential and temperature extracted from the measured particle multiplicity ratios. The results in Fig. 1 show the 
importance of the canonical treatment of strangeness conservation at SIS energies. Here, the canonical suppression 
factor can be even larger than an order of magnitude. For central Au-Au collisions at AGS or SPS energies this 
suppression is not relevant any more and the (GC)-formalism is adequate. In general, one expects that the statistical 
interpretation of particle production in heavy ion collisions requires the canonical treatment of strangeness conservation 
if the CMS collation energy y/s < 2 — 3 GeV. This is mainly because at these energies the f reeze-out temperature is 



still too low to maintain large-argument expansion of the Bessel functions in equation ( 2.12j ) 



At low temperatures one needs to take into account the width of resonances. This is because the number of pions 
coming from e.g. the decay of a A-resonance is increased by the width of the A. The approximation of the width by 
a delta function is therefore not justified since an appreciable number of particles come from the decay of resonances 



below the resonance mass. One therefore replaces the one-particle partition function in equations (|2.4|) and (2.5) by: 
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z r = n Xty t exp(6 fc ^B +qkVQ) [ dssK 2 (Vs/T)- t^tT 2 ^2 ( 2 - 15 ) 

where s m i n is chosen to be the threshold value for resonance decay and y/s ma x ~ m R + I\r. The normalization 
constant N is adjusted so that the integral over the Breit-Wigner factor gives one. 

Within the above model, the particle densities depend on four parameters: the chemical potentials, and hb, 
related with the (GC)-description of charge and baryon number conservation, the temperature T and the initial 
volume of the system appearing through the canonical treatment of strangeness conservation. Constraints on these 
variables arise from the isospin asymmetry measured by the baryon number divided by twice the charge, B/2Q. For 
an isospin symmetric system this ratio is simply 1, for Ni+Ni it is 1.04 while for Au+Au this ratio is 1.25. 

We are thus left with three independent parameters. For simplicity, the volume V will be identified with the volume 
of the system created initially in AA collisions estimated from the atomic number of colliding nuclei and from the 
impact parameter by using geometric arguments. In particular we use the relation of the volume parameter and the 
number of participating nucleons as it was indicated in Fig. 1. 

In the following section we will discuss to what extent the thermal model can be used to understand particle 
production in nucleus-nucleus collisions at SIS energies. 



III. COMPARISON WITH EXPERIMENTAL RESULTS 



This section is divided into three parts. In the first one we discuss the general trends found in thermal models. In 
particular we illustrate the sensitivity of particle ratios on the temperature, T, and on the baryon chemical potential, 
fiB- In the second part we discuss particle ratios measured in central collisions and compare experimental results 
obtained at SIS with the expectations of the thermal model. This avoids the problem that different particle species 
might originate from different impact parameter regimes as it is the case with inclusive studies. The impact parameter 
dependence is studied in the third part of this section. 



A. General Trends 



^From equations ( 2. 10| , 2. 12 2.13] ) one sees that at fixed temperature and chemical potentials the volume dependence 



of kaon multiplicity in the canonical and the grand canonical ensemble is as following: 

(N K )~{ V 2 . y ^ Q (c) _ (3.1) 

This effect is shown explicitly in Fig. 2 where the different particle ratios are calculated as a function of the radius of 
the system. In order to make the analysis more complete one needs, in addition, to take into account the contribution 
of resonance decays to particle production. We use all the known branching ratios as given in to calculate the 
particle multiplicities shown in Fig. 2. 



The ratio K + /ir + involves one strange particle, and according to equation (2.10) should show substantial dependents 
on volume. Indeed, it increases smoothly from zero up to the value given by the grand canonical ensemble. The 
i ncrea se is faster for higher temperature because the value of the argument x of the Bessel functions in equation 



(2.10) increases. In the K + /K~ ratio two strange particles are involved and the volume effect cancels out exactly. 
However, for very small values of the volume, there are a few non-strange resonances that decay into kaon pairs or into 
a kaon and a hyperon, this leads to the sharp rise in the value of this ratio for increasing volumes before it flattens off 
and becomes volume independent. The ratio <j)/i: + involves only non-strange particles and therefore is independent 
of the size of the system. The same behavior is expected for rj mesons which is treated as non-strange particle, too. 

Figure 3 evidences the locations of the freeze-out temperature Tf and of the freeze-out [i B yielding the various 
particle ratios. This figure shows also the sensitivity when varying their values within the limits occurring in the 
experimental results. Discussing these trends requires different treatment of strange and non-strange particles. 

The interesting feature of the deuteron to proton ratio, d/p, and the pion to proton ratio, tt + /p, is that they 
allow a good determination of the range of the thermal parameters. The it + /p curve in the (T — fis) plane shows 
a temperature saturation for large [1b which establishes the upper limit of the freeze-out temperature Tf. On the 
other hand, the d/p ratio fixes the range of the freeze-out n B as it shows a steep dependence on the temperature. 
In addition, as seen in Fig. 3, the variations of 20-30% on the value of the d/p ratio give a very similar range in the 
(T — hb) plane, making this observable particularly useful for the extraction of the freeze-out parameters. 
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In the SIS-energy range one expects, in the thermal model, a different dependence of the strange and non-strange 
particle yields on the volume of the system. Consequently, the ratio of strange to non-strange particle leads to a strong 
variation of thermal parameters with the system size. In Fig. 4 we show the K + /ir + ratio for different volumes. As 
expected, changing the volume implies a substantial modification of the curve in the (T — hb) plane corresponding 
to a fixed value of the K + /tt + yields. Thus, calculating the strange to non-strange particle ratio requires additional 
care of the system size. In our approach, however, the volume is not treated as an additional parameter but is rather 
related with the number of participating nucleons in AA collisions. For a given system size the K + /tt + ratio clearly 
determines the lower limit of the freeze-out temperature. The curves in Fig. 3 are calculated for R = 4 fm. The 
K + 1 K~ in Fig. 3 yield exhibits similar behavior as d/p, i.e. showing a very strong dependence on the temperature 
but a rather weak dependence on the baryon chemical potential, [Ib ■ This lines do not depend on the volume. 

An analysis of particle production in heavy ion collisions within a thermal fireball model requires two experimentally 
measured ratios to fix the freeze-out parameters, Tf, fi B , and knowledge of the number of participating nucleons, A part 
to establish the size of the fireball. The knowledge of more particle ratios allows to test the concept of a unique freeze- 
out time. In the following section we compare the predictions of the thermal model with experimental results for 
central A A collisions at SIS energies. 

B. Central collisions 

In central collisions, the number of participating nucleons is maximal. However, experimental results for zero impact 
parameter, 6 = 0, are not directly measured. However, for many experiments good-quality impact-parameter results 
are available and an extrapolation to b = can be performed. The results of this extrapolation are summarized in 
Table 1. We discuss below the different entries in this table. 

The results in Table 1 for pions are obtained from |l8| p(| . For the ratio ir + /p we used results from inclusive 
measurements since it was established that the pion multiplicity divided by the number of participants, M n /A par t, 
does not vary with A par t Jl8|-pT[. At Ni-Ni at 0.8 A-GeY n~ data are not available and n + = w~ /1.2 has been used 
to account for the isospin asymmetry. The K + results for Ni+Ni are from ]l4|] which are in very good agreement 
with p2| . The K + yield rises strongly with centrality as shown in |Q. At 1.0 and 1.8 A-GeY the K + /tt + ratio was 
obtained by extrapolating to b = 0. At 0.8 A-GeY we scaled with the inclusive K + /tt + ratio between 1.0 and 0.8 
A-GeY. The impact-parameter dependence for K~ and K + is nearly identical ]l3| and we used therefore the K + / K~ 
ratio of the inclusive measurements. The multiplicity of K + divided by the number of participants, M K + /A part , in 
Au+Au increases strongly with impact-parameter as shown in |l3|j23|1 . The ir~ /ir + ratios show - if at all - only a 
very slight increase with centrality J2(|^4| and Table 1 summarizes the results for inclusive studies. These values in 
f20y24f agree very well with the isobar model. The ratios of d/p obtained in Ni+Ni are from [^5[ and those for Au+Au 
have been taken from |26|. The ratios 77/71"° for Ni+Ni are from p7| , p8[ . At 1 A-GeY Ni+Ni has not been measured, 
yet Ar+Ca and Kr+Zr yield the equal ratios taking then also for Ni+Ni, at 0.8 A-GeV only Ar+Ca has been studied 
and this value is given in Table 1. The value at 1.8 A-GeY has been obtained by interpolation between 1.0 and 1.93 
A-GeY. For Au+Au the results from |27]] have been corrected for the increase with centrality according to p9[ . 

Figure 5 shows the lines in the (T — hb) plane corresponding to the measured particle ratios in Au-Au collisions 
at 1 A-GeV. The experimental errors are for simplicity not shown in the figure. All lines, except the one for 77/71"°, 
have a common crossing point around T ~ 50 and fis ~ 822 MeV. A value of R ~ 6.2 fm is needed to describe the 
measured K + /ir + ratio with the freeze-out parameters extracted for ir + /p, tt + /tt~ and d/p. This radius corresponds 
to Ap ar t ~ 330 participating nucleons and is compatible with the one expected for central Au-Au collision. 

The strangeness suppression due to the canonical treatment of the conservation laws is very clear in the comparison 
of the thermal model with the Au-Au 1 A-GeY data. Using the grand canonical formulation of the strangeness 
conservation one would get the value K + /tt + ~ 0.04 which overestimates the data by more than an order of magnitude. 
This shows that the conditions for thermal particle phase space at SIS energies are far from the grand canonical limit. 

In Fig. 6 we show the corresponding results in the (T — fis) plane for Ni-Ni collisions at three different incident 
energies: 1.8, 1.0 and 0.8 A-GeV. As for the Au-Au data we see in Fig. 5 that all particle ratios, besides 77/71-° can 
be described by the same values of the freeze-out parameters. We notice the smaller radius of ~ 4 fm which is 
compatible with the smaller size of Ni. The measured K + /K~ ratio leads to a band in the (T — hb) plane which 
barely misses the common crossing point for the mean values of 7r + /p, d/p and K + /ir + . However, taking into account 
the experimental errors on the above particle ratios leads to a common, narrow (T — /is)-band which contains also 
the line corresponding to the upper experimental limit for K + /K~ ~ 30. Allowing for a drop in the K~ mass as 
proposed in model calculations pp|-|32]| leads to a shift of the K + /K~ band towards the left in Fig. 6, thus leading 
to a better agreement with the other data. New experimental results on the K + /K~ ratio also for other collision 
systems are needed to clarify this open question. 
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The values for the of d/p ratio shown in Fig. 6 for Ni-Ni collisions at 0.8 A-GeV were obtained by extrapolating 
the experimental measurements at E/A=1.06, 1.45, 1.93 summarized in Table 1 using a linear and a polynomial fit 
giving the value d/p = 0.4, 0.43. As one can see the 10% deviation on d/p ratio does not modify substantially the 
freeze-out line in the (T — ps) plane. 

The results of Figs. 5, 6 show that for central Au-Au and Ni-Ni collisions the particle ratios, n + /p, K + /ir + , ir + /ir~ , 
K + 1 K~ and d/p, lead to a common crossing point in the (T — /ib) plane. The appearance of the common freeze-out 
for all these particles is strong support for chemical equilibrium of these particles in the thermal model. 

In Fig. 7 we show the results for the freeze-out temperature and chemical potential corresponding to different 
incident energies. It can be seen that the freeze-out temperature Tf increases with E/A whereas pr B shows the 
opposite trend. In the energy range considered this dependence can be approximated by a straight line. The freeze- 
out parameters are also seen in Fig. 7 to be different in Ni-Ni and Au-Au collisions even for the same incident energy. 
It is interesting to remark that at a given incident energy the extracted values of Tf are smaller for central Au-Au 
collisions than for central Ni-Ni ones. 

Using the freeze-out parameters relevant for SIS energies we can compare the results with the previous findings for 
AGS (3^]34|] and SPS p5|-|38[ energies. The chemical freeze-out points in the (T — ps) plane for relativistic heavy ion 
collisions are shown in Fig. 8 together with the predictions for e + e~ and p-p collision s 139(1 . Connecting these points 
leads to a unique freeze-out curve in the (T — ps) plane [j40|j4l| l. It has been shown [|l2| that the energy per hadron 
along the freeze-out curve (before the final decay of the hadron) is approximately 1 GeV. 

It is a common feature of Figs. 5 and 6 that the freeze-out line for rj production does not cross the common 
chemical freeze out extracted for all the other particle species. This deviation does not arise from the selection of 
central collisions. Only for Au-Au collisions we have chosen the 77/77° ratio of central collisions. In the Ni-Ni system - 
due to the lack of experimental results - inclusive values are given. For central Ni-Ni collisions the discrepancy would 
therefore increase. The results in 29 show that the ij yield rises more than linearly with A part . This disagrees also 
with the expectations from the thermal concept discussed along with Fig. 2. 

The problem of rj abundance in the thermal model could have several origins. One possibility would be due to 
the sequential freeze-out for different particle species. Here, 77 could be produced earlier with higher temperature, 
roughly corresponding to the value obtained from the transverse momentum slope parameter of 80 MeV. In this 
case, as seen in Fig. 6, the large value of the 77/77° ratio can be understood. In recent works j2^,^,Q experimental 
evidence is given that high-energy pions are emitted earlier than those with lower energies. This points towards a 
span in freeze-out times but does not affect the particle ratios which are dominated by low-energy pions. The concept 
of sequential freeze-out, however, might lead to a problem with the interpretation of all other particle ratios which 
are otherwise well described by the model. We stress that the same thermal model when applied to AGS and SPS 
energies explains particle productions there with a single set of freeze-out parameters for all particles including 77 at 
SPS energies. It seems rather difficult to argue the appearance of a sequential freeze-out at SIS and its absence at 
AGS and SPS energies. 

The discrepancy of the thermal model on the level of r\ production could possibly be related to the hidden strangeness 
content of the r\ meson. The crucial difference in the thermal model interpretation of particle production at SPS, 
AGS and SIS energies is due to the canonical strangeness suppression. In our formalism the 77 is considered as a 
non-strange particle. It is conceivable that due to hidden strangeness there should be corrections to the rj yield from 
the canonical strangeness conservation. We do not know, yet, however, how these corrections could be included in a 
consistent way. The same arguments would apply also for the production of </>. Yet, <j)/K~ (inclusive) seems to fit 
into the common crossing shown for Ni-Ni at 1.8 A- GeV. The large rj yield could well be of dynamical origin which 
apriori cannot be explained by the thermal freeze-out model applied in this work. Thus, we leave the discrepancy of 
rj production as an unsolved puzzle in the thermal model. 

Concluding this part, we find that the thermal model gives a consistent description of both strange and non-strange 
particle production in central nucleus - nucleus collisions at SIS energies with the exception of the rj/ir° ratio. There 
is a serious problem with the understanding of this ratio in terms of the thermal model. Theory and experiment 
deviate by a factor of 2 to 8, depending on the incident energy. The same problem has been indicated recently in the 
context of dynamical model for pion and rj production |16| . We would like to point out that recent transport- model 
calculations reproduce the measured 77 and 77° spectra lft5| . We would also like to draw attention to the systematics 
introduced by V. Metag in which pions and 77's exhibit a common trend while K + , K~ and <ft show lower yields Km . 

To check further the consistency of the model with the experimental data the impact parameter dependence of 
particle production and particle spectra are discussed in the next section. 
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C. Impact Parameter Dependence 



As already mentioned, the multiplicity of K + divided by the number of participants, K + / A part , increases strongly 
with the centrality while the corresponding ratio for pions, K + /A part , is constant [ fblf . Consequently, the pion yield 
is proportional to the number of nucleons in the initially created fireball while the multiplicity of K + scale with A par t 
like K + ~ A part with a > 1. The experimental results on K + /ir + and K + /A part ratios in Au-Au collisions from Ref. 
p3f are shown in Fig. 9. The canonical treatment of strangeness conservation predicts the yield of strange particles 
to increase quadratically with the number of participants (see eq. (2.14)). An additional complication is due to the 
possible variation of the freeze-out parameters with centrality. 

In Fig. 9 the dashed-line describes the results of the thermal model under the assumption that both Tf and n B 
are A part independent. One sees that already under this simple approximation the agreement of the model and the 
experimental data is very satisfactory. The small deviations between the model and the data can be accounted for by 
the variation of the freeze-out parameters with A part . 

We have calculated the possible change of Tf and [i B with A part from the experimental data on K + /tt + and d/p 
ratios measured for two different values of A part (see Fig. 10). The resulting freeze-out parameters are shown in 
Fig. 11. One can see that the variation of Tf and [J B with A part are small. It is interesting to note that peripheral 
collisions yield higher temperatures than central collisions. The results from central Ni — Ni collisions fit perfectly 
into the trend with A part . 

In Fig. 9 the dashed-dotted lines were obtained parameterizing the Tf and \J B dependence on A part using the small 
variation in freeze-out parameters obtained from Fig. 11. The shape of the experimental data is now well reproduced. 

The rather small variation of the freeze-out temperature with impact parameter shown in Fig. 11 comes as a 
surprise since the experimental results on the apparent inverse slope parameter T app of particle yields shows a strong 
dependence on A par t . The freeze-out temperature of ~ 50 MeV derived from our analysis is also substantially lower 
than the value previously obtained from the particle spectra |l8| - ^0| , p5| , ^4 47 4qj . 

We now turn our attention to the differential cross sections for particle production. In Fig. 12 we show the variation 
of the inverse slope parameters observed for various particle species |l§| , |25]| . Yet, adopting our result, that there exists a 
common, impact-parameter independent, freeze-out temperature for all particles one needs to show that the variation 
in slope results from radial flow varying with impact parameter. We extract the values of the corresponding flow 
parameters (3 using the Siemens- Rasmussen formula p9[ |. The resulting values for the flow velocity, (3, are summarized 
in Fig. 12. With the freeze-out temperature Tf ~ 50 MeV the spectra of p, K + and 7r + are, within the experimental 
uncertainty, well explained with the same (3 as shown in Fig. 12. The variation of (3 with impact parameter turns out 
to be very large. 

The A part dependence of various particle yields (Y ~ Ap art ) has been discussed in ]r| and the exponent a has been 
suggested to be related with the difference in total energy needed for the production of the studied particle and the 
one available in TV TV collisions. In the presented concept the exponent a is only due to volume-dependent strangeness 
suppression. In this picture the exponent for K + and for K~ production is predicted to be equal while in the frame of 
the energy argument jl8| a should be higher for K~ than for K + due to different production thresholds (neglecting 
in- medium mass modifications). It is clear that any non-isotropic emission pattern is beyond the scope of this model. 
However, - as naively expected - the slopes of K + and of K~ do not have to be equal as different resonance decays 
contribute. 

The results of this section show that the thermal model gives a consistent picture for particle production in AA 
collisions at SIS energies. In addition to the correct predictions for particle yields the model can also explain the 
Apart dependence of K + cross section. Assuming radial flow as an origin of the shape of particle spectra one can 
understand quantitatively the inverse slopes with a common temperature and radial flow velocity for all particles. 
The chemical and thermal freeze-out seem to be very close at SIS energies. This is because the particle multiplicities 
and momentum spectra can be explained with the same temperature. 



IV. SUMMARY 



In summary, the thermal model provides a remarkably consistent description of the experimental data in the 
GSI/SIS energy range. The abundances of particles, K + , K~ ,p, d,ir + and tt~ (with the notable exception of 77's) 
seem to come from a common hot source with a surprisingly well defined temperature, T « 50, 54, 70 MeV, and 
baryon chemical potential [is ~ 825,805,750 MeV for central Ni-Ni at 0.8, 1.0, 1.8 A-GeV and correspondingly 
T w 52 MeV and [is ~ 822 MeV for central Au-Au collisions at 1.0 A-GeV, as can be seen clearly from Figs. 5,7 and 
8. These temperatures are lower than the ones observed in the particle spectra but here again a common explanation 
is possible in terms of hydrodynamic flow. Flow differentiates between heavy and light particles since they acquire 
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the same boost in velocity but of course very different kicks in momenta. This is clearly seen in the GSI/SIS data 
shown in Fig. 12 which also summarize the transverse flow velocity for pions, kaons and protons. 

The common crossing points exhibited in Figs. 5, 6 are a very strong argument for chemical equilibrium. The 
deviations found for the production of r\ mesons clearly ask for an explanation. 
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TABLE I. 
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VI. FIGURE CAPTION 

Fig.l Canonical strangeness suppression factor (see text). 

Fig. 2. Dependence of various particle ratios on the radius parameter R. 

Fig. 3 General trends for ir + /p, d/p, K + /n + and K + /K~ freeze-out curve calculated with R = 4.0 fin and isospin 
asymmetry corresponding to A/2Z = 1.04. 

Fig. 4 Dependence of K + /tt + freeze-out line with the radius R .. 

Fig. 5 Freeze-out lines corresponding to different particle ratios measured in Au+Au collisions at 1 A- GeV. 
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Fig. 6 T versus /ib for central Ni+Ni collisions from 0.8 to 1.8 A-GeV. 

Fig. 7 Freeze-out Tf and yJ B as a function of incident energy as extracted from the common crossing in Figs. 5, 6 - 
neglecting the results from n/ir . 

Fig. 8 Freeze-out parameters Tf and \x B for SPS, AGS and SIS energies. 

Fig. 9 Measured K + /ir + and K + multiplicity per A part as a function of A part for Au+Au at 1 A-GeV together with 
two calculations (see text). 

Fig. 10 Freeze-out parameters T and fi B for different A part . 

Fig. 11 Variation of freeze- out parameters T and [Ib as a function of A part . 

Fig. 12 Apparent temperatures for p,n + , K + (data from Jjj|/j and calculated flow velocities as a function of A par t 
for Au+Au at 1 A-GeV. 
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